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Let us suppose that customers arrive at a counter in accordance with a 
Poisson process of density X. The customers are served by a single server in 
order of arrival. The service times are identically distributed, mutually in- 
dependent, positive random variables with distribution function II (x). 
Suppose that after being served each customer either immediately joins the 
queue again with probability p or departs permanently with probability q 
(V + Q = I)- In this paper we shall determine for a stationary process the 
distribution of the queue size as well as the Laplace-Stieltjes transform and 
the first two moments of the distribution function of the total time spent in 
the system by a customer. 

I. INTRODUCTION 

Although the problems discussed in this paper arose in the theory of 
telephone traffic, we use the terminology of queues. Thus instead of calls 
and holding times we shall speak about customers and service times 
respectively. 

Let us suppose that in the time interval (0,<x>) customers arrive at a 
counter in accordance with a Poisson process of density X. Denote by 
r n (n = \ , 2, ■ ■ ■ ) the arrival time of the nth customer. Then the inter- 
arrival times t„ +1 — t„ (n = 0, 1, • ■ ■ ; r = 0) are identically distrib- 
uted, mutually independent random variables with distribution function 

F(r) = l l - e ~ XX iix = °> (1) 

1{X) \0 if u; < 0. K) 

The customers are served by a single server in order of arrival. The 
server is idle if and only if there is no customer in the system. The service 
times are supposed to be identically distributed, mutually independent, 
positive random variables with distribution function H(x), and inde- 
pendent of the input process. Suppose that after being served each cus- 
tomer either immediately joins the queue again with probability p or 
goes away permanently with probability q where p + q = 1 . The event 
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that a customer returns is independent of any other event involved and, 
in particular, independent of the number of his previous returns. 

The process denned above is said to be of type [F(x),H(x),p]. If p = 0, 
then there is no feedback. 

Let us denote by £(£) the queue size at time t, that is, £(Z) is the num- 
ber of customers (either waiting or being served) in the system at time t. 
Let £„ denote the queue size immediately before the arrival of the nth 
customer, that is, £„ = £(t„ — 0). Denote by 0„ the total time spent in 
the system by the nth customer. 

We are interested in finding the distribution of %(t) for a stationary 
process, i.e., when %(t) has the same distribution for all t ^ and the 
distribution of 0„ for a stationary process, i.e., when 8 n has the same 
distribution for every n = 1,2, • • • . 

It is easy to prove that the limiting distribution lim P{£(0 = k\ 

t-*oo 

(k = 0, 1, • • • ) exists and is independent of the initial state if and only 
if $(<) has a stationary distribution and the limiting distribution is 
identical with the stationary distribution. Similarly the limiting dis- 
tribution lim P{ d„ ^ x\ exists and is independent of the initial state if 

and only if 6„ has a stationary distribution and the limiting distribution 
is identical with the stationary distribution. 

Throughout this paper we use the following notation : 

*(s) = r e-"dH(x) (SR(s) ^ 0) (2) 

for the Laplace-Stieltjes transform of H(x), 

ot r = f K 'x'dH(x) (r = 0, 1, •••) (3) 

for the rth moment of II (x), and 

a = [" xdH(x) (4) 

for the average service time, i.e., a = «i . 

Further denote by H*(x) the distribution function of the total service 
time of a customer. We have 

H*(x) = gJlp^Huix) (5) 

t=i 

where H k (x) denotes the kth. iterated convolution of H(x) with itself. 
For qp k ~ l (k = 1, 2, ■ • • ) is the probability that a customer joins the 
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queue k times, and if he joins k times, his total service time is equal to the 
sum of A- mutually independent random variables each of which has the 
distribution function H(x). If we introduce the notation 

**(«) = f" <T tx dH*(x) (SR(s) = 0), (6) 

then by (5) we obtain that 

•n.)-, gj^iwr-j^j. (7) 

Let 

a r * = [* x T dH*(x) (r = 0,1, ••■) (8) 

•'o 
and a* = on*. By (7), a* = a/q, 

^-S + fe", (9) 

q ([- 

and in general a* can be obtained by the following recurrence formula 
*r* = a r + ?t(-)^_y* (10) 

II. THE STATIONARY DISTRIBUTION OF THE QUEUE SIZE 

If we know the stationary distribution of the queue size for a process 
of type [F(x),H(x)0], then that for a process of type [F(x),H(x),p] can 
be obtained immediately. 

Theorem 1: IfXa < q, then the process [l-(t),0 = t < =o} has a unique 
stationary distribution P{£(0 = j\ = Pj*(j — 0, 1, • ■ •) and for | 3 | = 1 

// Xa ^ 5, J/ien a stationary distribution, does not exist. 

Proof: To find the distribution of the queue size we may assume with- 
out loss of generality that the customers join the queue only once and 
are served in one stretch ; however, their service time is equal to the total 
service time that they would have if they were served in the original 
manner. Accordingly the distribution of the queue size for the process 
of type [F(x),H(x),p] is the same as for the process of type 

[F(x),H*(x),0]. 
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For the latter a stationary distribution {Pj*} exists if and only if Xa* < 1, 
that is, Xa < q, and is given by A. Y. Khintchine's formula: 

UH») = I> V = ° - X "* )(1 ~^ ,(X(1 - Z)) . (12) 



i=o 



¥*(X(1 - z)) - z 



(Cf., e.g., Ref. 1 or 2.) This proves (11). 

Remark 1. Denote by B*(r = 0, 1, ■ • • ) the rth binomial moment of 
{Pi*}, that is, 

B T * = 1LfyPj*- (13) 

If a r+ i is finite, then B r * is also finite. We have B * =1, 

d * x2 " 2 * _li * - A 2 a 2 + 2X«i(l - Xai) n .v 

Bl = 2(1 - Xa:*) + Xai 2( 9 - Xa.) (14) 



and for r = 2, 3, 

X r 

r! |_£i (1 - Xa 

where 



. „ X ^ XV! r r fi XM F I 

£i (1 — Xa*)" X —i (1 — Xa*)" J 



(15) 



y = y nla,*''c»*'' ••■«,+>«'■ (16) 

"' W^+ik-. jilj.1- ■ - j.l(2l)"(3!)" • • • ((n + DO'" 
ii+2ia+' ■■+nj„=n 

The proof of (14) and (15) can be found in Ref. 3. 

III. THE STATIONARY PROCESS 

Let us denote by x« the time needed to complete the current service 
(if any) at the instant t = t„ - 0, i.e., immediately before the arrival of 
the nth customer. If £„ = 0, then x» = 0. It is easy to see that the vector 
sequence {£„ , x» J n = 1, 2, • • •} is a Markov sequence. We shall prove 
that if Xa < q, then {£„ , Xn] has a unique stationary distribution, 
whereas if Xa ^ 5, then a stationary distribution does not exist. For a 
stationary sequence {£„ , x») introduce the following notation: 

Pi - Fie. -i) 0" = o, 1, ■••), (17) 

Pi(s) = P{xn ^ -r, £„ = 3] U = 1, 2, • • • ) (18) 



and 



M = f C-dPAx) (j = 1,2, •■•)• (19) 

•'n 



y 
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Theorem 2: If\a < q, then the Markov sequence {£„ , x« ; n = 1, 2, • • •} 
has a unique stationary distribution, for which P, = P * defined by (11) 
and 



U(s,z) = £ n } (s)z' 



j=i 



0-t) 



(20) 
Xg(l -jOj^a) -*(A(1 -z))] 

[z - (q + p*)*(X(l - *))][« - X(l - *)] ' 



// \<x ^ 5, JJken {£„ , Xn} /'«s no stationary distribution. 

Proof. First consider the process of type [F(x),H*(x),0]. It is proved 
in Ref. 3 that in this case {£„ , Xn} has a unique stationary distribution 
if and only if \a < q. Namely P{£„ = j\ = P,* (j = 0,1, • • • ) given by 
(11) and the generating function corresponding to (20) is given by 

_,„, (1 - Xa)te( l - z)[**(s) -**(X(1 - z))] r91 v 

U *M= [«-»*(X(l -z))][s -A(l-z)] • (21) 

The distribution of the queue size is the same for both the types 
[F(x),H(x),p] and [F(x),H*(x),0). 

The only difference between the process of type [F (x) ,H (x) ,p] and 
[F(x),H*(x),0] is that in the latter x» the remaining length of the cur- 
rent service at the arrival of the nth customer is replaced by the remain- 
ing part of the total service time of the customer just being served at 
the arrival of the nth customer. The time added to Xn is independent of 
the queue size and has the distribution function 

H(x) = qt,pH k (x), (22) 

k=0 

because the probability that a departing customer will join the queue 
k more times is qp k and in this case the additional total service time has 
the distribution function H k (x); H (x) = 1 if x ^ and H (x) = if 
x < 0. The Laplace-Stieltjes transform of (22) is 

$(,) = r e~" dm*) = q E p k Ms)] k = : ^ . (23) 

Accordingly we have 

U*(s,z) = U(s,z)Ms), (24) 

whence (20) follows. 

Remark 2. If %{t) denotes the time needed to complete the current 
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service (if any) at time t(x(t) = if £(0 = 0), then the vector process 
(£(£), x (0;0 ^ t < oo } is a Markov process. {£(t)iX(0} has a stationary 
distribution if and only if Xa < q and it agrees with the stationary 
distribution of { £„ , Xn\ • 

IV. THE STATIONARY DISTRIBUTION OF 0„ 

If the joint distribution of £„ and x» is known, then the distribution of 
0„ is determined uniquely. If J £„ , x»l has a stationary distribution, then 
every 6 n (n = 1,2, •■ ■) has the same distribution. In case of a stationary 
process let 

E{ e - s6n ) = $(s) (5R(s) ^ 0). (25) 

Theorem 3: If Xa < q, then 0„ has a unique stationary distribution 
P{0 n S %}, which is given by the following Laplace- Stieltjes transform 

Hs) = qT,p k - 1 U k (s,l) («R(8) ^ 0), (26) 

where 

C/i(s,«) = Po^(s + X(l - «)) 

+ U(s + X(l - z),(q + P*)*(* + X(l - z))) (27) 

/or 9i(s) ^ and | 2 | ^ 1, Po = 1 — Xa/o, J7(s,2) is defined by 
(20), and 

t/ fc+ i(s,2) = *(s + X(l - z))U k (s,(q + P*)*(s + X(l - z))) (28) 

/or fc = 1, 2, ■ • • . 

Proof: The probability that a customer joins the queue exactly k 
times (including the original arrival) is qp _1 (fc = 1,2, • • •)• Denote 
by d n {k) the total time spent in the system by the nth customer until his 
fcth departure (if he joins the queue at least k times). Denote by £ n w 
the queue size immediately after the fcth departure of the nth customer. 
Let 

U k (s,z) = E{exp [se n {k) ]z uik) \ (k = 1, 2, • • •)• (29) 

We can easily see that for a stationary sequence {£„ , x>«j 
Ui(s,z) = P<&(s + X(l - z)) 

+ E lUs + X(l -«))(« + P»)1*(« + X(l - «))]' 

3=1 ( 30 ) 

= Po*(s + X(l - z)) 

+ C/(s + X(l - z), (q + pz)*(s + X(l - z))) 
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where P = 1 — \a/q and U(s,z) is denned by (20). Now we shall 
prove that (28) holds for k = 1, 2, • • ■ . Under the conditions {„ (k) = j, 
0„ = .r and that after the fcth service the roth customer joins the queue 
again, the difference 0„ ( * +1) - 0„ ai is equal to the length of j + 1 ser- 
vices, the distribution function of which is// j+i(x). If 0„ a+1> — 9 n {k) = y, 
then f„ " is equal to the sum of two independent random random vari- 
ables : the first is the number of customers arriving at the counter during 
the time interval of length y, which has a Poisson distribution with 
parameter \y. and the second is the number of returning customers, which 
has a Bernoulli distribution with parameters,/ and p. Thus 

E|exp [-.s0/+V" a ' +n I r» W - j, On™ = x) = e-(q + V z)>. (31) 
/ exp[-sy - X(l - z)y] dH j+l (y) 

= «-"(« + j*)'[*(i + X(l -z))\ J+ \ 
or 

E{exp[-«fl. (H ' , V* < * +l) |r. tM ,0 

(32) 
■ *( S + X (l - z) ) exp [-aO.Kfl + P*)*(s + X(l - j))/"*** 

and unconditionally 

17* +X («,*) - *(s + X(l - z))U k (s,(q + p*)*(s + X(l - 2))) (33) 

which proves (28). It is to be noted that stationarity has been used only 
in the determination of Ui(s,z). The recurrence relation (28) is valid for 
any process. Finally, 

Ejexp [sd n ]} = q f) p^Ejexp [-s0 B (t) ]} = q f] p M ^( S) l) (34) 

which was to be proved. 

V. THE MOMENTS OF 0„ 

Although it seems very complicated to find a closed formula for $($), 
the moments of 0„ can be determined explicitly. We shall prove that 
$>(s) = 4?(s,l ) where *(s,z) satisfies a functional equation. This observa- 
tion makes it possible to find explicit formulas for the moments of 0„ . 

Theorem 4- If Xa < q and 0„ has a stationary distribution, then 

b(«,i = xa, + wi-xa,) (35) 

2{q - Xai) 
provided that a-> is finite, ami 
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EIO - s q " 2q 



6(g - Xai) 2 [g 2 - g(2 + Xa,) + X«i] 

• {2g[6X«i 3 - 6ai 2 - 6Xaia 2 + 3a 2 + Xa 3 ] (36) 

- [12X«! 3 - 12ai 2 - 6Xaia 2 + 2X 2 aia 3 - 3XW]}, 

provided that a 3 is /ira'te. 
Proof: Let 

*(«,*) = <ziV _1 ^-(s,z) (37) 

fc=i 

where E/*(s,2) is denned by (27) and (28). By using the recurrence 
relation (28) we obtain that 

*(«,«) = qUi(s,z) 



(38) 
+ p*(s + x(i - z))*(«,(g + ps)*(s + x(i - z))). 



Let 



_ /a' w «(«,«)\ . (39) 

J \ ds'ds' / s =o, 2 =i 

If we form $„■ by (38) for * + i =" r (* = 0, 1, • • • , r), then we obtain 
r + 1 linear equations for the determination of *,-,• . These equations can 
be solved successively for r = 1, 2, • • ■ . By (26) 

E{0 n r ) = (-1)'** (r = 0, 1, •••) (40) 

for a stationary process. 

VI. A PARTICULAR CASE 

If we suppose, in particular, that the service times have an exponen- 
tial distribution 

*(,)-(!-'"" J"* J (41) 

v ' (0 for a; < 0, 

then H*(x) = 1 - e~ mx for * ^ 0, ¥(«) = m/(m + «) and ¥*(a) = 
M9/(M2 + s). In this case (27) and (28) reduce to 

"*<*) - ^./xo-.) p > ('• ; /.Vx/T - J (42) 

f or k = 0,1, • • • , where 

^-O-sVO-s)- (43) 
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By (42) we obtain that 



U k (s,l) = 



i) = V ml 
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(44) 



Ok(s) - b k (s) 

where a k (s) and b k (s) (k = 0, 1, • • • ) are given by the following matrix 
equation 

M + X + s 



fik(s) 



h(s) 



X 



(45) 



Now by (26) the Laplace-Stieltjes transform of P{0„ ^ a;} is 



*(•) - 1 -- 



(i-A)£ 



qp 



Zi a k (s) — b k (s) ' 



In this particular case 



and 



Effc) - 



E[0„ 2 } = 



I 



»q - \ 



2 M (2g - r/) 



(m</ - X) 2 [ M (2 9 - 5 2 ) -X(l - q)Y 



(46) 



(47) 



(48) 



VII. ACKNOWLEDGMENTS 

I am indebted to Dr. H. O. Pollak and Mr. J. P. Runyon, who called 
my attention to the problem discussed in this paper. I would like also to 
thank Dr. W. S. Brown for programming an IBM 7090 computer at the 
Bell Telephone Laboratories, Murray Hill, to determine the second 
moment of 0„ which is given by formula (36). It is most remarkable that 
Dr. Brown's method made it possible to use an electronic computer for 
performing algebraic operations. 

APPENDIX — BY W. S. BROWN 

Calculation of the Second Moment of 0„ 

If Xa < q and if 8„ (the total time spent in the system by the nth 
customer) has a stationary distribution, then the moments of 6„ are 
determined by (38)-(40) as explained in Section V. The first moment 
can be calculated by hand without serious difficulty. The calculation of 



514 THE BELL, SYSTEM TECHNICAL JOURNAL, MARCH 1903 

the second moment, with the aid of an IBM 7090 compute]-, 4 is de- 
scribed in this appendix. 

It is convenient to replace z by I = 1 - z and Ui by W = qUi/ 
(q — X«i). Then the rth moment of 6„ is 

f} r =(-lY*ro (49) 

where 

The function *(s,0 is implicitly defined by the equation 

*(«,<) = (g - \ ai )W(s,t) + pHs + X0*(*,«(M)) (51) 

whore* 



with 



., v ^ (-l) r a r s r 

u(s,t) = 1 - (1 - pt)+(s + \t) 
W(s,t) = ^(s + \t) + £(s + X<,X«(*,0)7'(m(8,0 



a; ?/ 

71/ \ _ Ml - to) 

KW> 1 - « - (1 - pw)^(Xw) ' 



(52) 



(53) 



This last pair of equations can be rewritten in the more useful form 

X(l - w) 



(54) 



r(«) = - 



g — X(l — pu)<p(\u) 



where 



r+l r+l r 

_ .r — y _ ^ * r-fc 



C r (.r,*/) = ^— = 2-X y 

X — 7/ t -0 

, v 1 - ^(:c) A (-l) r a, + iX r 

**>-— r-"S (r+l)! 



(55) 



* For convenience we have assumed that all of the service moments a r [see (3)1 
arc finite. However for the calculation of /3 r it is clearly sufficient to require only 
the finiteness of a r +i. 
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It is now clear that 

*(0) = <* = 1 
5(0,0) = -ax 

T(Q) = - (56) 



q — Xai 



W(0,0) = 



q — Xai ' 



so from (49)-(51) 



/3 = *oo = $(0,0) = 1 (57; 



as is required by the definition of the zeroih moment. 

Now suppose all of the quantities $ij for i -f- j < r, where r is some 
positive integer, have been calculated and are expressed as rational 
functions of X and p (or q) and the service moments <n- • Then by differ- 
entiation of (51 ) we can obtain a system of r + 1 linear equations in the 
r ■+■ 1 unknowns, <l>, 7 with i + j = r. These equations will also contain 
the quantities <!>,•_,■ with i + j < r, which can be replaced by their known 
values. The solutions of this linear system will again be rational func- 
tions of X and p (or q) and the service moments a k . Theoretically this 
procedure permits the calculation of arbitrarily many of the moments, 
but in practice the calculations are extremely lengthy. The reader may 
wish to try the first moment as an exercise. 

We shall now outline the computation of the first two moments, ft 
and fa . The first step was to fix the time scale by setting X equal to 
one. The next step was to express $(s,<) as a Taylor series with coeffi- 
cients fyj and expand (51) to second order in s and t. Since rational 
function operations were not yet available, all denominators had to be 
eliminated by suitable premultiplications. The result obtained by the 
computer was a gigantic polynomial in s, t, <f> o , $01 , $10 , $02 , $11 , $20 , 
ai , as , a,i , and q. We chose to view it as a polynomial in s and t. Setting 
$00 = 1, all the terms independent of both s and t vanished as expected. 

The equations for $01 and $10 , obtained by setting the coefficients of 
t and s to zero, are 

2A<b Ql + 5 = 

(58) 
2C*io + 2D*oi + E = 
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where 

A = (q- ai )[q - g(«i + 2) + orj 

B = -{2ai - 2 ai - ao)[q* - ?(<*i + 2) + «i] 
C = -q(q - ai) (59) 

D = —ai(q - l)(q - en) 
E = (2«! 2 - 2«i - a 2 )f(«i + l)g - «i]. 
These expressions were factored by hand. From (58) we have 

B 



(60) 



_ BD - AE 
*" 2IC - ' 

so, using (59) 

2ai — 2ai — <*2 / fi1 \ 

$01 = $10 = 7w S • ^ l) 

2(q - a x ) 

Thus the mean time spent in the system (waiting time plus service 
time) is 

fc _ -* M = «, (Lz«) + «• . (62) 

\q — ax} 2{q — on) 

In the absence of feedback (p = 0, q = 1 ) this reduces to 

rf"-«+*rHj (63) 

where the first term is the mean service time and the second term is 
the familiar expression 5 for the mean waiting time in a single server 
system. 

The equations for $02 , $11 , and <J> 20 , obtained by setting the coeffi- 
cients of t 2 , st, and s 2 to zero (and replacing $10 by $01), are 

A 13 3>02 + Bi*oi + Ci = 

^22$ll + A 23 <J> 02 + #2*01 + C? = (64) 

A 81*20 + A 32 $ll + ^33*02 + -S 3 $01 + C 3 = 
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where 

-4,3 = 6(g - ai)V - q\2ay + 3) + q( ai 2 + 4a, + 3) - ai ( ai + 2)] 

^ 22 = -12(g - a,) 2 [ 5 2 - g(a, + 2) + oj 

A n = — 12ai(? - 1)(? - ai) 2 (g - a, - 1) 

2 (65) 

-43i = Qq(q — on) 



A 32 = 12a,(g - l)(g - a,) 2 

-433 = 6a, 2 (g - l)(g - a,) 2 

and 

£, = 6(g - l)(g - a,) 2 [4tfa, - (2a, 2 + 4a, + a 2 )] 

B 2 = 12(2 ~ 1)(2 - ai) 2 [2gai - (2a, 2 + 3a, + a 2 )] (66) 

#3 = -6(g - l)(g - a,) 2 (2a, 2 + 2a, + a 2 ) 

and finally 

& = - 2? 4 (6a, 3 - 6a, 2 - 6a,a 2 + 3a 2 + a,) + 3 ? 3 (8a, 4 -8a, 2 a 2 

— 8a, — 6a,a 2 + 2a,a 3 — a 2 2 + 6a 2 + 2a 3 ) — 2<? 2 (6a, 5 
+ 12a, — 6ai a 2 — 6ai — 21a, 2 a 2 + 3a, 2 a 3 — 12a, 2 

— 3a,a 2 2 + 7a,a 3 - 3a 2 2 + 9a 2 + 3a 3 ) + g(12a, 5 (67) 

— 18a, a 2 + 2a,' a 3 — 3a, a 2 — 12a, a 2 + 10a, 2 a 3 

— 12a," — 9a,a 2 " + 12a,a 2 + 10a,a 3 — 6a2 2 ) 

+ a,(6a,"a2 — 2a, a 3 + 3a,a 2 " — 6a,a 2 — 4a,a 3 + 3a 2 2 ) 
Ci = 2g'(12a, — 18a, — 12a, a 2 + 6a, 2 + 6a,a 2 + 2a,a 3 + 3a 2 
+ a») - ? 2 (24a, 5 + 12a, 4 - 24a, 3 a 2 - 48a, 3 - 48a, 2 a 2 
+ 8a, 2 a 3 + 12a, 2 - 6a,a 2 2 + 30a,a 2 + 12a,a 3 + 3a 2 2 
+ 18a 2 + 6a 3 ) + tf(24a, 5 _ I2 ttl 4 - 36a,a 2 + 4a, 3 a 3 (68) 

— 12a, — 6a, 2 a 2 — 12a, a 2 + 14a, a 3 — 9a,a 2 2 + 30a,a 2 
+ 12a,a 3 — 3a 2 ") + ai(12a, 2 a 2 — 4a, 2 a 3 + 6a,a 2 2 

— 12a,a 2 — 6a,a 3 + 3a 2 ") 



(69) 
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C 3 = -q{\2ax - 12cn 4 - Ua^ - 6«i 2 a 2 + 2 ai 2 a 3 + 6ai« 2 

+ 2aia» + 3a 2 2 + 6« 2 + 2a 3 ) + gai(12ai* - 12a/ 

- 18ai a 2 + 2ctiat - Zona? + 4aia 3 + 12a 2 + 4a 3 ) 

4-ai 2 (6aia-> — 2«iar 3 + 3a 2 — 6a 2 — 2a 3 ). 

These expressions were factored by the computer using the "divide if 
divisible" subroutine and a long list of trial divisors including all factors 
appearing in (59). The complete solution of (64) was obtained in three 
passes through the computer with some assistance from the author 
between passes. On the first pass the quantity 

Bi*oi + Ci 

was calculated [using (61) for * ra ] and an attempt was made to divide 
it by each of the factors of A 13 . The division by the cubic factor was 
successful, and there followed the result 

* M= J^G (70) 

Q{q - ai)- 

where 

F = 6ai 3 - 6ai 2 - 6ai« 2 + 3a 2 + a ? , 

G = 12aj 3 - 12ai 2 - 6«ia 2 + 2ai« 3 - 3a 2 \ 

On the second pass the quantity 

i4 23 $ 2 + #2*01 + Ci 

was calculated [using (70) for $02 and (61) for <*><»] and an attempt was 
made to divide it by each of the factors of A 22 and by the numerator of 
$02 • The latter division was successful, and the result 

[2qF - G][q - g(«i + 3) + gj (y9) 

11 12 (q - ai) 2 fa* - ff(ai + 2) + aj 

was thereby obtained. On the third pass the quantity 

[q - q(oti + 2) + aJ[iliAi + ^33*02 + # 3 <i>oi + C 3 ] 

was calculated (the first factor being introduced in order to cancel the 
corresponding factor in the denominator of *n) and an attempt was 
made to divide by each of the factors of A Zi . Only the division by q 
was successful. This yielded the final result 

_ (2qF - G)(g* - 2q) , . 

ft " * 2 ° - 6 ( g _ „ lW _ g(ai + 2) + «d U * ; 



A SINGLE-SERVER QUEUE WITH FEEDBACK 519 

where F and G are the expressions defined in (71). Note that $02 appears 
as a factor in <t>n and <J> 2 o . In the absence of feedback (p = 0, q = 1) 
(73) reduces to 

„<«> = ^_ + »' + ^L_^ (74) 

1 - «! 2(1 - aj- 3(1 — ai) 

which is the correct result. To see this we use the law of composition 
to obtain 

/3'2 = a-> + 2aia>i + coo (75) 

where cti and « 2 are the service moments while a>i and co 2 are the waiting 
moments. The latter are given in Ref. 5 as 

_ "- 
Ul " 2(1 - «i) 

(76) 

402 ~ 2(1 - ax) 2 3(1 - ax)' 

Substituting (7G) into (75) we obtain (74) as expected. To obtain 
(30), substitute (71) into (73), replace each a t . by \ h a k , and divide the 
resulting expression by X 2 . 
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